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Abstract 

Using data collected near the D*^D^ peak production energy E cm = 4170 MeV by the CLEO-c 
detector, we search for Df exclusive hadronic decays involving u. We find B{Df — > tt + uj) = 
(0.21 ± 0.09 ± 0.01)%, B{D+ -► Tr+Tr ^) = (2.78 ± 0.65 ± 0.25)%, B{D+ -► rt+^-w) = (1.58 ± 
0.45 ± 0.09)%, B(D+ -► 7r+r/Lj) = (0.85 ± 0.54 ± 0.06)%, B{Df -► K+cu) < 0.24%, £(L>+ -► 
K+vr ^) < 0.82%, B(L>+ -► i^+vr+vr-aj) < 0.54%, and B(D+ -► if+r/u;) < 0.79%. The upper 
limits are at 90% confidence level. 
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The inclusive u> yield, Df — > ooX, is substantial (6.1 ± 1.4)% This is very surprising, 
as the only exclusive decay mode involving u that has been observed is — > 7r + u;, 
with a branching fraction of B(Df — > h + uj) = (0.25 ± 0.09)% [2J. There is lots of room for 
more exclusive hadronic decays involving uo. The study of uo production in decays 
is of interest in shedding light on mechanisms of weak decay and their interplay with long- 
distance (nonperturbative) physics [3|. Here we present a search for several exclusive 
hadronic decays involving u. In particular, we consider final states with one, two, and three 
pions: — > — * 7r + 7r°u;, and — > 7r + 7r + 7r~c<j. We also search for — > 7i + r]LJ, 
as this has been suggested jl] as possibly being a large decay mode. Finally, we search 
for modes in which one of the 7r + from the above-mentioned decays is replaced by a K + : 
— > K + u, — > K + 7t°uj, — > K + 7i + 7i^uj, and — > K + rjuj. These last four modes 
would be Cabibbo-suppressed, and hence are not expected to be large. 

In this study we use 586 pb -1 of data produced in e + e~ collisions at the Cornell Electron 
Storage Ring (CESR) and collected by the CLEO-c detector near the center-of-mass (CM) 
energy a/s = 4170 MeV. At this energy the cross-section for the channel of interest, D* + D~ 
or D+D*~, is approximately 1 nb [4]. We select events in which the D* decays to D s ^ 
(94% branching fraction [2}). Other charm production totals ~7 nb and the underlying 
light-quark "continuum" is about 12 nb. 

The CLEO-c detector is a general-purpose solenoidal detector, which is described in detail 
elsewhere [5], @, 0, S]. The charged particle tracking system covers a solid angle of 93% of 
47r and consists of a small-radius, six-layer, low-mass, stereo wire drift chamber, concentric 
with, and surrounded by, a 47-layer cylindrical central drift chamber. The chambers operate 
in a 1.0 T magnetic field. The root-mean-square (rms) momentum resolution achieved with 
the tracking system is approximately 0.6% at p = 1 GeV/c for tracks that traverse all layers 
of the drift chamber. Photons are detected in an electromagnetic calorimeter consisting of 
about 7800 CsI(Tl) crystals [oj. The calorimeter attains an rms photon energy resolution of 
2.2% at E 1 = 1 GeV and 5% at 100 MeV. The solid angle coverage for neutral particles in 
the CLEO-c detector is 93% of 4ti. We utilize two particle identification (PID) devices to 
separate charged kaons from pions: the central drift chamber, which provides measurements 
of ionization energy loss (dE/dx), and, surrounding this drift chamber, a cylindrical ring- 
imaging Cherenkov (RICH) detector, whose active solid angle is 80% of An. The combined 
PID system has a pion or kaon efficiency > 85% and a probability of pions faking kaons (or 
vice versa) < 5% |9j. The detector response is modeled with a detailed GEANT-based llOj 
Monte Carlo (MC) simulation, with initial particle trajectories generated by EvtGen [ll| 
and final state radiation produced by PHOTOS [13]. The modeling of initial-state radiation 
is based on cross sections for D* Df production at lower energies obtained from the CLEO-c 
energy scan [i| near the CM energy where we collected the sample. 

Here we employ a double-tagging technique, the same as the technique that is used in the 
D+ inclusive decay analysis [l|. Single-tag (ST) events are selected by fully reconstructing 
a D~ , which we call a tag, in one of the following three two-body hadronic decay modes: 
D~ — > K° S K~ , D~ — > <f)7t~ and D~ — > K*°K~ . (Mention of a specific mode implies the 
use of the charge conjugate mode as well throughout this Letter.) Details on the tagging 
selection procedure are given in Refs. flHQ. The tagged D s candidate can be either the 



primary D s or the secondary D s from the decay D* — > ^D s . We require the resonance 
decays to satisfy the following mass windows around the nominal masses Q|: K$ — > 7r + 7r~ 
(±12 MeV), (p -> K + K~ (±10 MeV) and K*° -> K + n~ (±75 MeV). The momenta of all 
charged particles utilized in tags are required to be 100 MeV/c or greater to suppress the 
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FIG. 1: The mass difference AM(D S ) = M(D S ) — m£> g distributions in each tag mode. We fit the 
AM{D S ) distribution (points) to the sum (solid curve) of signal (double-Gaussian) plus background 
(second degree polynomial, dashed curve) functions. 



slow pion background from D*D* and D*D decays (through D* — > tcD). 

The reconstructed invariant mass of the D s candidate, M(D S ), and the mass recoiling 

against the D s candidate, M Tecoi i(D s ) = <J (E — E Ds ) 2 — (p — PdJ 2 , are used as our pri- 
mary kinematic variables to select a D s candidate. Here (E , p ) is the net four-momentum 
of the e + e~ beams, taking the finite beam crossing angle into account, po a is the momentum 
of the D s candidate, Ed s = \l m< D s + Pd s , an d mr> s is the known D s mass [2]. We require the 
recoil mass to be within 55 MeV of the D* mass [2j. This loose window allows both primary 
and secondary D s tags to be selected. We also require a photon consistent with coming 
from D* — > decay, by looking at the mass recoiling against the D s candidate plus 7 
system, M Teco ^{D s ^) = ^ (E — E Da — E^) 2 — (p — Pd s — P^) 2 - For correct combinations, 
this recoil mass peaks at m^ s , regardless of whether the candidate is due to a primary or a 
secondary D s . We require |M recoil (D s 7) — m^J < 30 MeV. 

The distributions of AM(D S ) = M(D S ) — rri£> a) the invariant mass difference of D s tag 
candidates, after applying the two recoil mass requirements for each tag mode are shown 
Fig. [TJ We use the tag invariant mass sidebands to estimate the backgrounds in our signal 
yields from the wrong tag combinations (incorrect combinations that, by chance, lie within 
the AM(D S ) signal region). The signal region is \ AM(D S )\ < 20 MeV, while the sideband 
region is 35 MeV < \AM(D S )\ < 55 MeV. To find the sideband scaling factor, the AM(D S ) 
distributions are fit to the sum of double- Gaussian signal plus second-degree polynomial 
background functions. We have 18586 ± 163 ST events that we use for further analysis. 

In each event where a tag is identified, we search for our signal modes recoiling against 
the tag. Charged tracks utilized in signal candidates are required to satisfy criteria based 
on the track fit quality, have momenta above 50 MeV/c, and angles 9 with respect to the 
beam line, satisfying |cos6*| < 0.93. They must also be consistent with coming from the 
interaction point in three dimensions. Pion and kaon candidates are required to have dE/dx 
measurements within three standard deviations (3<t) of the expected value. For tracks with 
momenta greater than 700 MeV/c, RICH information, if available, is combined with dE/dx. 

We identify 7r° candidates via 7r° — > 77, detecting the photons in the Csl calorimeter. 
To avoid having both photons in a region of poorer energy resolution, we require that at 
least one of the photons be in the "good barrel" region, |cos# 7 | < 0.8. We require that 
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the calorimeter clusters have a measured energy above 30 MeV, have a lateral distribution 
consistent with that from photons, and not be matched to any charged track. The invariant 
mass of the photon pair is required to be within 3cr (a ~ 6 MeV) of the known 7r° mass. 
The 77 candidates are formed using a similar procedure as for 7r° except that o ~ 12 MeV. 
The 7T° and 77 mass constraints are imposed when 7r° or 77 candidates are used in further 
reconstruction. We reconstruct u candidates in the u — > 7r + 7r~7r° decay mode. 

Mode-dependent requirements on numbers of charged kaons and pions are applied on the 
signal side. For example, we require there must be exactly one charged kaon and two charged 
pions for the D£ — > K + uj,uj — > 7r + 7r~7r° mode. No extra tracks are allowed on the signal 
side. The best n° (or 77) candidate is selected based on the pull mass (number of standard 
deviations of measured mass from true mass). For — > n + r]uj, a veto on 77' — > 7r + 7r~77 
has been applied to remove the dominant background contribution from the — > 7i + 7r°rj' 
decay. 

The double-tag (DT) yields are extracted from the 7r + 7r~7r° invariant mass distribution 
after requiring that both the tagging D s and signal D s invariant masses be in the D s nominal 
mass region (20 MeV mass window on the tag side and 30 MeV mass window on the signal 
side due to it or 77 on the signal side). The uo mass signal region is \M ir + ir - ir o —mj\ < 20 MeV, 
while the sideband region is 40 MeV < \M n +^-^o — m w | < 80 MeV, where M^+^-^o is the 
7T + 7r~7r° invariant mass and m w is the nominal mass of uj |2J. 

The invariant mass distributions of u candidates for the first four modes are shown in 
Fig. [2] and the yields are given in Table [H For — > tt + uj, — > 7t + tt°uj, and — > 
7t + 7t + 7t~uj, clear signals are found in the data. Note also the peaks from 77 — > 7r + 7r~7r° and 
— > 7T + 7r _ 7r , corresponding to known D s decays to 77 plus pions and plus pions. 

For .D+ — > 7i + i]uj, we observe 7 events in the signal region, and 5 events in a sideband 
region twice as wide. From Binomial statistics, the probability of observing 7 or more events 
in the signal region and 5 or fewer events in the sideband region, out of a total of 12 events, 
if there is no true signal, is 6.6%. Thus we have "evidence" for — ► tt + i]uj, but cannot 
claim "observation". We quote both a value for the branching fraction and an upper limit 
on it. 

We find no significant evidence for any of the modes with kaons, and therefore set upper 
limits on their branching fractions. The numbers of events from u mass signal and sideband 
regions for these four modes are given in Table [B Upper limits are calculated using Poisson 
statistics, and allowing for uncertainty in number of background events. 

We determine the efficiency for detecting our tags (cst) using Monte Carlo samples in 
which one D s decays into the tag mode, and the other decays generically. The single-tag 
efficiency weighted over the three tags is (26.63 ± 0.14)%. We determine the efficiency for 
detecting both tag and signal (6dt) using Monte Carlo samples in which one D s decays 
into a tag mode, and the other D s decays into a signal mode. For the three-body decay 
modes, we assume three-body phase space. For — > tt + 7t°uj, we also consider — > 
p + u,p + -> tt+tt . For D+ -> K + n u, we also consider D+ -> K* + u,K* + -> K + tt°. For 
— > 7t + 7t + 7t~uj, we assume 7t + 7t + 7t~uj phase space and estimate a systematic error by 
using _D+ — > 7i + p°u,p° — > 7r + 7r _ . Double-tag detection efficiencies are given in Tabled! 

The absolute branching fractions (£> mo dc) are obtained from the ST yield (A^st) and DT 
yield (A^dt) without needing to know the integrated luminosity or the produced number of 
D^DJ pairs, 

iV DT e S T / 1 \ 

Smodc = X . (1) 
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FIG. 2: Invariant mass distributions of u candidates for the first four modes: — > it + uj, — > 
7T + 7r°w, — > 7T + 7r + 7r _ w, and — > 7r + r]io. The solid lines (red online) indicate the to mass 
signal region and the dashed lines (blue online) indicate the u mass sideband regions. Peaks from 
7] — > 7T + 7T~7r and (j) ~^ 7r + 7r~7r° are also evident. 



The .D+ — > 7r + 7r°CG> decay might come from — > p + (7r + 7r°)u;. We fit the corresponding 
7r + 7T° invariant mass distribution to the sum of phase space ir + n°u MC and p + uj MC. The fit 
result suggests that 0.52 ± 0.30 of the — > 7i + 7i°u) decay come from the — > p + u, p + — > 
7i + 7T° decay, as shown in Fig. [31 

We have considered several sources of systematic uncertainty. The uncertainty associated 
with the efficiency for finding a track is 0.3%; an additional 0.6% systematic uncertainty for 
each kaon track is added [^]. The relative systematic uncertainties for tt° and 77 efficiencies 
are 4.0%. Uncertainties in the charged pion and kaon identification efficiencies are 0.3% 
per pion and 0.3% per kaon For the — > tt + 7t°uj mode, the relative contribution of 
Df — > p + uj, p + — > 7r + 7T° is determined to be 0.52 ± 0.30 from the fit. We use the central 
value of that ratio to calculate the efficiency and take the error as a systematic uncertainty. 
All Monte Carlo efficiencies have been corrected to include several known small differences 
between data and Monte Carlo simulation. Upper limits have been increased to allow for 
the systematic errors. 

The branching fractions and upper limits are listed in Table [TTJ In summary, we report 
first observations of — > tt + 7t uj and £)+ — > 7t + 7t + tt~u) decays. The branching fractions are 
substantial. We find evidence for the — > n + r]u) decay. Our measurement of — > -k + uj 
decay is in good agreement with the PDG value |2j, and of comparable accuracy. The 
sum of branching fractions of these four observed modes is (5.4±1.0)%, which accounts for 
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TABLE I: Observed yields. Here Ag g is the observed event number from u mass signal region, 
iVsd is the scaled event number from u mass sideband regions, and Ag s is the sideband-subtracted 
signal yield. Double-tag detection efficiencies (cdt) are listed in the last column. The efficiencies 
include sub- mode branching fractions [2|, and have been corrected to include several known small 
differences between data and Monte Carlo simulation. Errors shown are statistical errors only. 
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FIG. 3: 7T + 7r° invariant mass distribution. We fit to the sum of phase space 7r + 7r°u; MC and 
p + u MC. The fit result suggests that 0.52 ± 0.30 of the Df — » tt + tt°uj decay come from the 
— > p + uj,p + — > 7r + 7r decay. The points are the data and the superimposed line is the fit. 

TABLE II: Branching fractions and upper limits. Uncertainties are statistical and systematic, 
respectively. 
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most of the D s inclusive u decays (6.1±1.4)% [lj]. We also report the first upper limits on 
D+ -> K+ou, D+ -> K+ttV D+ -> K+n+n-u, and L>+ -> K+r/w decays. 
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